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Abstract 

We discuss a maximally localized Wannier function approach for constructing lattice models 
from first-principles electronic structure calculations, where the effective Coulomb interactions 
are calculated in the constrained random-phase-approximation. The method is applied to the 3d 
transition metals and a perovskite (SrVOa). We also optimize the Wannier functions by unitary 
transformation so that U is maximized. Such Wannier functions unexpectedly turned out to be 
very close to the maximally localized ones. 
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I. INTRODUCTION 



In a class of materials often referred to as correlated materials, the electronic structure 
is characterized by a set of partially filled narrow bands across the Fermi level. Seen from 
the atomic site, one has a set of partially filled shell of localized orbitals typically of 3d 
or 4f character. Many of the electronic properties of the material are determined by the 
correlations among the localized electrons living in the partially filled band or shell. It is 
therefore physically well motivated to map the original complicated many-electron problem 
to a model consisting of the localized orbitals and a few additional orbitals. By eliminating 
the high-energy states ("downfolding") the long-range bare Coulomb interaction is screened 
to a short-range interaction at low energy. Since the screened interaction is short range, only 
on-site interaction or the Hubbard U is often taken into consideration in the model. This 
is the physical idea behind the well-known Hubbard model or Andersen impurity model. 
One would then wish to have a set of well localized orbitals or Wannier orbitals that span 
the same Hilbert space as that of the states that form the narrow bands. In this way the 
Hubbard U will have small off-site matrix elements, which may be neglected. Practical 



procedures to construct the mode 



s starting from first-principles calculations have been a 



subject of interest for a long time l|, 3, sl, 

In this work, we focus on Wannier orbitals using the method developed by Souza, Marzari 
and Vanderbilt based on the minimization of the quadratic extent of the orbitals. An 
alternative, equally promising approach is to use the Wannier orbitals of Andersen jrj]. While 
the former is a "post-processing" method, i.e., the Wannier orbitals are constructed after 
generating the Bloch wave functions, the latter may be termed "pre-processing" method 
because the Wannier orbitals are constructed before diagonalization of the Hamiltonian that 
yields the ba.d structure g In th. sense, the .atte. sche,„e „,ay be „,o.e advantageous 
than the former. On the other hand, the former is more general because it does not depend 
on any particular band-structure method. Comparison between the two Wannier functions 
for some selected materials can be found in Ref.Q]. 

Apart from the use of Wannier orbitals in constructing lattice models, there are many 
other applications. In particular, a close connection with Berry's phase [s, ^ has stimulated 



intensive works recently 



12 



13l ? ]. Related work can be found not only in condensed 



matter physics but also in chemistry, where the concept of localized molecular orbitals is 
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very useful for understanding chemical bonding as well as for visualization. The idea of 
constructing localized molecular orbitals goes back to the early sixties by the maximization 
of the Coulomb energy of the molecular orbitals ij] or the minimization of the quadratic 



extent of the molecular orbitals 



15| . The general problem of transforming a set of Bloch 



states to a set of well localized orbitals is therefore one of important methodological problems 
in condensed matter physics and chemistry. 

Another important issue in the downfolding procedure is how to determine effective in- 



teraction parameters. A widely used method is constrained LDA (cLDA) 1, ^16| and a 
recently proposed scheme based on the maximally localized Wannier function 17|] may be 



useful for applications to complicated structures. On the other hand, cLDA is known to 
yield unreasonably large values of U in some cases (e.g. late transition metals). This arises 
from technical difficulty in including part of the self-screening processes between localized 
electrons leading in some cases to a larger value of U 18|. Another method for estimating 
effective interaction is the random phase approximation (RPA). We can find early trials 
along this line in Ref. 19|, l20t]. Later on the constrained RPA (cRPA) scheme was invented 
jsl. The cRPA method has several advantages over currently available methods. It allows 
for a precise elimination of screening channels, which are to be included in a model Hamilto- 
nian, without modifying the one-particle dispersion of the model. In addition, the effective 
screened interaction as a function of r and r' can be calculated independent of the basis 
functions. We will use this method in the present work. We can also find other proposals 
in literature such as a hybrid method between cLDA and cRPA 2l(| and linear response 
approach [22 1. 

Our long-term goal is to construct a first-principles scheme for calculating the electronic 
structure of correlated materials. As is well known, the local density approximation (LDA) 



23| in density functional theory (DFT) [24] often has difficulties when applied to such 



systems. Attempts to improve the LDA have resulted, among others, in the LDA+U la. 
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26(, the LDA+DMFT (Dynamical Mean-Field Theory [271 l28l |) |29| and more recently in 



the newly developed GW+DMFT method 30|, ImI]- In Ref. 31| it is shown how the Hubbard 



U for real materials can be determined self-consistently within the scheme. In these methods 
it is crucial to have well localized orbitals representing the Hilbert space of the partially filled 
correlated bands since the screened Coulomb interaction U is usually assumed to be purely 
on-site. This is especially the case in the DMFT method, where the lattice problem is 
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mapped to an impurity problem 32 1. 

The purpose of the present work is to demonstrate the usefulness of the Wannier orbitals 
and the feasibility of performing many-body calculations with a unified Wannier basis, in- 
dependent of the starting band structure. To this end we have calculated the Hubbard U 
within the cRPA scheme using the maximally localized Wannier basis and compared the 
results with independent calculations j^, [3] in the linear muffin-tin orbital (LMTO) basis 



331]. The reasonably good agreement between the two sets of results gives us confidence as 



to the usefulness of the scheme. 

As mentioned earlier, it is highly desirable to construct a set of Wannier orbitals that 
minimize the off-site Coulomb interaction or equiyalently maximizes the on-site U. We fol- 
low the method of Edmiston and Ruedenberg ij], which was proposed for molecules, and 
developed a practical procedure to maximize the U parameter for periodic crystals through 
unitary transformation in real space. Application to transition metals shows that the effect 
of maximization is tiny if we start the optimization from the maximally localized Wannier 
functions. 



II. METHOD 

The Wannier function with band index n at cell R is defined by 

M = / e^^'-^'l^ilV'^ , (1) 

where is the associated Bloch function which can be expanded as a linear combination 

of the eigenfunctions of a mean-field Hamiltonian as 

l^i^) = E^-(k)|^™k). (2) 

m 

In practical implementations, Kohn-Sham wavefunctions may be used for IV'mk)- In the 
maximally localized Wannier function scheme , the coefficients Wmn(k)'s are determined 
such that the quadratic extent of wavefunctions 

^ = ^{{Vno\r'^Wno) - \{^no\AVno)?) , (3) 
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TABLE I: Energy window 1 that limits states to be included in constructing Wannier functions, 
and window 2 to specify localized orbitals {ipd} for cRPA. The "2nd band" for the window 2 means 
the second lowest one in the bands with strong 4s and 3d character. Energies are measured from 
the Fermi level. 





Window 1 


Window 2 


Sc 


[-3.0 eV, 5.0 eV] 


[2nd band, 4.05 eV] 


Ti 


[-4.0 eV, 5.0 eV] 


[2nd band, 3.85 eV] 


V 


[-4.0 eV, 5.0 eV] 


[2nd band, 4.15 eV] 


Cr 


[-5.0 eV, 4.0 eV] 


[2nd band, 2.85 eV] 


Mn 


[-5.0 eV, 4.0 eV] 


[2nd band, 1.50 eV] 


Fe 


[-5.0 eV, 4.0 eV] 


[2nd band, 1.20 eV] 


Co 


[-5.0 eV, 3.0 eV] 


[2nd band, 0.55 eV] 


Ni 


[-7.0 eV, 3.0 eV] 


[2nd band, 0.25 eV] 


SrVOs 


[-10.0 eV, 5.0 eV] 


3 t2g states 



is minimized. For this purpose, we introduce an energy window and optimize Wmn(k) with 
limiting m to the states inside the window. The parameters for this window ("window 1") 
are Usted in Table [B The Wannier function is more localized as the energy window is larger, 
since optimization is done in wider Hilbert space. We found, however, that the (screened) 
Coulomb interaction is not sensitive to the choice of the energy window unless the window 
is too wide. 

The idea of cRPA is to define an effective interaction Wr by excluding screening processes 
that are included in an effective low-energy model with Wr as the effective interaction 
(Hubbard U). To this end we divide the Hilbert space into two parts: localized states that 
form the projected bands {ipd} and the rest {ipr}- The polarizability P is then divided into 
two as P = Pd + Pri where Pd includes transitions between {il)d} only and P^. is the rest 
of the polarization. The effective interaction in the reduced space Wr is defined so that 
Wr[l — PdWr\~^ yields the RPA fully screened interaction W = v\l — PvY^, where v is the 
bare Coulomb interaction. It can be shown that such Wr is given by 

Wr{uj) = V[l - Pr{uj)v]-^ , (4) 
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where spatial coordinates are omitted for simplicity js^. This Wr, after multiplying some 
localized functions and integrating over space, can be interpreted as the frequency-dependent 
Hubbard U. We refer the static value of Vrr(O) as the Hubbard U used in model Hamiltonians. 

In the following sections, we compute matrix elements of Wr in the Wannier basis. The 
calculation starts with a conventional LDA electronic structure obtained by the full-potential 
LMTO method. The polarizability is computed using Kohn-Sham eigenf unctions and eigen- 
values. It can be computed efficiently for an arbitrary number of frequencies at essent ially the 
cost of just one frequency [3J]. Other technical details are found elsewhere |35|, [Sg, l37|, |38 |. 



We use a 8 X 8 X 8 k mesh for transition metals and a 4 x 4 x 4 k mesh for SrVOs for the Bril- 
louin zone integration. Optimization of Wannier functions is done following the procedure 



in Ref . 



l6|. 



III. CRPA WITH THE MAXIMALLY LOCALIZED WANNIER FUNCTION 

Let us start with paramagnetic nickel as an example. We first construct five Wannier 
orbitals having strong 3d character. Once Wmn(k) is determined on a k mesh, maximally 
localized Wannier functions are obtained by Fourier transform, from which the Hamiltonian 
Hjnn(R) = {'^mo\H\ipnii) IS rcduccd as well. By Fourier transforming i/m„(R) back to k 
space and diagonalizing it Q], we can project out narrow bands (FigJT]). 

The next step is to compute the screened Coulomb interaction Wr{r,r']uj) in cRPA and 
take the matrix elements in the maximally localized Wannier basis: 



Wr{ni,n2,n3,n4;R;u) = j y (/7;^o(r)^„,o(r)lK(r, r'; ^)(/?;^i^(rV„4R(r')rfWr' . (5) 

At this point it is worth pointing out that the effective screened interaction Wr{r,r';uj) 
calculated using the cRPA method is completely independent of the choice of basis functions. 
The matrix elements are of course dependent on the choice of the orbitals (/9„o(r)- The on- 
site diagonal elements, Wr{n,n,n,n,Il = 0; uj) {n = I,-- - ,5), are split by crystal field 
effect, however the splitting is negligibly small. Figure [2] shows the average of the five, 
where fully screened Coulomb interaction is also shown for comparison. Wr is close to the 
bare Coulomb value (dot-dashed line) at high energy, since screening effect is minor. As 
the frequency decreases down to around 30 eV, the screening becomes effective and Wr 
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FIG. 1: Bandstructure of paramagnetic nickel in LDA (solid thin lines) and projected bands (dotted 
thick lines). 

decreases rapidly. At lower energy Wr is weakly energy dependent again and reaches 2.8 eV 
at u; = 0. These features are the same as the previous calculations [3|, [l8|, though the values 
are slightly smaller in the present result. The difference may be ascribed to the difference 
in spacial extent of the orbitals ipnoiT")- In the previous calculations, the effective screened 
interaction Wr{r, r'; u) is calculated within the LMTO-ASA scheme and the orbitals v5„o(r) 
are taken to be the truncated partial waves, i.e., the heads of the LMTO or the solutions 
of the Schrodinger equation inside the atomic sphere, rather than the LMTO basis. They 
are normalized and completely confined to the atomic spheres. In the present calculations, 
Wr{r, r'; u) is calculated using the full-potential LMTO (FP-LMTO) scheme but the orbitals 
(fno are the maximally localized Wannier functions, which have tail extending outside the 
central cell. These orbitals are thus more delocalized and consequently the matrix elements 
of Wr are smaller. The difference in Wr{r, r'; uj) arising from the difference between LMTO- 
ASA and FP-LMTO is probably less significant. 

The above informations, projected band structure and on-site values of Wr, are key 
ingredients for constructing effective models. However, it is not sensible to construct a 
Hubbard model by simply adding the former as the kinetic term and the static value of 
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FIG. 2: On-site screened Coulomb interaction of paramagnetic nickel in the maximally localized 
Wannier basis as a function of frequency. Fully-screened interaction (closed circles), result from 
cRPA (open circles), and bare Coulomb interaction (dotted line). 

the latter as the interaction term. In fact, excitation spectra for that Hamiltonian do not 
reproduce the solution of the original Hamiltonian. This is because (i) the kinetic term is 
renormalized during the downfolding process, and (ii) Wr is energy dependent and long- 
ranged. However, there is an approximate way to construct a Hubbard Hamiltonian with a 
static interaction ^|. 

Figure[3]^a) shows the diagonal element of the on-site Wr in the static limit (U) for a series 
of transition metals. Comparing the results with those obtained from previous calculations 
presented in (b), the trend is the same: As the atomic number increases, the U increases in 
the early transition metals, while it decreases in late transition metals. We also computed 
U in SrVOs (Table [IT|) . In this system, there are three t2g states near the Fermi level and 
they are isolated from other bands. Thus, there is no ambiguity for dividing the space into 
{i^d} and {i^r}- The value of U is computed to be 3.0 eV, which is again smaller than the 
previously calculated value of 3.5 eV. We emphasize again that in the previous calculations 
the effective screened interaction Wr is calculated within the LMTO-ASA scheme and the 
choice of the orbitals in calculating the matrix elements of Wr are truncated partial waves. 
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FIG. 3: On-site screened Coulomb interaction U for transition metals from cRPA (open circles) and 
fully screened U from RPA (closed circles), (a) Results when the matrix elements ofU = Wr{uj = 0) 
are taken in the maximally localized Wannier function basis, and (b) in the truncated partial waves 
or the heads of the LMTO-ASA basis. 

which are confined within the atomic sphere and thus more locahzed compared with the 
Wannier orbitals used in the present calculations, leading to a larger value of U. The U 
value obtained for SrVOs should be similar for CaVOa, LaTiOs, and YTiOs perovskites. 

The off-diagonal (exchange) elements are also important quantities. In Fig|4](a), 
{Wr{n,m,m,n,Il = 0; uj))n^rn in Ni is shown as a function of frequency, where average 
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TABLE II: On-site Coulomb (U), exchange (J), and off-site Coulomb ([/') energy in SrVOa 
obtained by cRPA. 



U 



3.0 eV 



U' 



0.45 eV 



J 



0.43 eV 



is taken over n and m. In contrast with the Coulomb term, the exchange term is weakly 
energy dependent and does not show significant change at around 30 eV. We may under- 
stand this behavior as follows. / (i^r'Wr(r, r'; a;)y9*gj^(r')(y9„^R,(r') is a screened potential of 
a charge density p(r') = V2^3R(r')v^n4R(r')- accordance with known observation, the po- 
tential arising from an exchange charge density (713 7^ 77,4) is not well screened because it 
has no mono pole (zero spherical average) in contrast to the case of = n^. At the onset 
of the plasmon excitation at around 30 eV, a perturbing charge is highly screened and the 
screening is electron-gas-like, which is highly effective for a mono pole. At lower energy (5-6 
eV), however, atomic-like screening in the form of 3d-3d transitions as well as 3d-4p takes 
place, which is relatively effective in screening a multipole charge distribution, resulting in a 
significant decrease of J. As can be seen in the case of Fe and Ni, J is reduced considerably 
when 3d-3d screening arising from is included whereas in Cu, where there are essentially 
no 3d-3d transitions, since the 3d band is fully occupied, the main screening channels come 
from 3d-4p transitions and J varies less strongly compared to those of Ni and Fe for the 
fully screened case. Figure [5] shows the static value of the exchange term, J. We find that 
J does not depend on the element significantly and its value is around 0.5 eV. 

Another important information is non-locality of the interactions. The bare Coulomb 
interaction v as a function of R = |R| is long ranged (Figj6](a)), and v between the nearest 
neighbor cells is about 1/4 of the on-site value. On the other hand, the screened interaction 
shows much faster damping and the value at the nearest neighbor {U') is 0.1 eV, which 
is much smaller than the on-site value of U=2.7 eV. The values of U' for other transition 
metals are shown by open circles in FigjSl 

The approximately parabolic variation on both U and U' across the series may be quali- 
tatively understood in terms of 3d band filling. The largest polarization inside the 3d band 
(Pd) corresponds roughly to half- filling. Since this is eliminated when calculating U and U', 
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FIG. 5: On-site screened exchange interaction J (triangles) and off-site screened Coulomb inter- 
action U' (open circles). 



ully screened interac- 
since the screening 



they peak around the middle of the series. This is on contrast to the 
tion W which is almost a constant across the series. As discussed in 
is metallic it does not depend much on the element: there are always enough electrons to 
screen a perturbing charge. Moreover, since the screened interaction is rather localized, it 
is not sensitive to the extent of the orbitals used in taking the matrix elements in Eq.(l5]). 

The U' is larger in some elements. For example, U' of Cr is as large as 1.0 eV, which 
would not be negligible. In such cases, one would wish to reconstruct Wannier function 
so that off-site interaction is as small as possible. This possibility is discussed in the next 
section, where a procedure to maximize the on-site U is derived and applied. 



IV. MAXIMIZING THE ON-SITE U PARAMETER 



A. Formulation 

We follow closely the method in Ref. [l^ and apply it to the case of periodic crystals. We 
use the convention that repeated indices are summed. Let us define a unitary transformation 
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FIG. 6: Diagonal elements of Coulomb interaction of Ni as a function of i? = |R| (in units of 
lattice constant), (a) bare Coulomb interaction, and (b) Screened Coulomb interactions in cRPA 
(open circles) and in RPA (closed circles). 



XnR = ^nR + S^nK 



(6) 
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From the unitarity of T one has, to first order 



(7) 



S^nR — '/^n'R'T'n'R'.nR, Tn'R',nR + "^n'R'.nR — . (8) 

Consider a change in U to first order in Stpnii- 



U = J2ix:nXnR\U\x:nXnR) (9) 

nR 

= (¥'n'R''y^nR|t^bnR¥'nR) <'R',nR + {¥>*nRVn'R'\U\ipl^ipnR) Tn'R'.nR 
+ ('y^nR<y^nR|C^kn'R'<y^nR) + {^*nB,^nR\U\(pltiVn'R') T„'R',„R 

= 2 {^n'-R'VnIi\U\(plji(PnIl) <'R',nR + 2 (V'nRV'n'R' 1 V^hRV^hr) Tn'R',nR • (10) 

This provides an expression for the change oi U as a function of independent parameters 
"^n'R'.nR '^n'R'.nR- Using t*/r/ „r + t„r,^„/r,/ = One Can rewrite this expression as follows. 



SU = -2 {v*n>n>VnR\U\(p*„ji(p„R) rnR,„/R/ + 2 (<yf*R'yfn'R'|t^bnR<yfnR) T„'R',nR 

= 2 [{v*nR^n'R'\U\(pl^ipnR) " (V^nR^^n'R' | ^^n'R' '/'r^'R' )] ^n'R',nR 

^ '^FnR,n'R'^n'R',nR , (H) 

where we have defined an anti-Hermitian matrix (F'^ — —F) 

KR,n'R' = {VnRfn'R'\U\(p*^^fnR) " {^nR^n'R'\U\ip*^,^,^n'R') = F*'R',nR ■ (12) 

We now choose 



Tn'R'.nR — ^-^n'R',nR , (13) 

which ensures that 



5U{e) = 2£F+,^„,j,,F„,R, „R ^ . 
14 



(14) 



The procedure is then the following. Construct the matrix 

T = e"^ , (15) 

which is unitary because F is anti-Hermitian. To calculate T we diagonalize F with eigen- 
vectors {nK\a) and eigenvalues fa'- 

Tn'wM^) = {n'I{!\a) e'f- {a\nR) . (16) 
One now obtains a new basis and calculates a new [/ as a function of e. 



U{e) = (XnRXnR|f/|XnRXnR) 

= {'^iR.iVjB.2\U\(pl^.^(piIl^) ^*Ri,„R^jR2,nRrfcR3,nR^«R4,nR • (17) 

One varies e until U{e) reaches a maximum and chooses that new basis set XnR. that max- 
imizes U{e). The procedure is then repeated until convergence is achieved. In practice, we 
solve for the eigenvectors and eigenvalues of iF, which is Hermitian. The eigenvalues of F 
are then given by —i times the eigenvalues of iF. The starting orbitals v9„r, are chosen to 
be the maximally localized Wannier orbitals but other choices are also possible. 

The key quantity in the present formulation is the anti-Hermitian matrix F, which is 
defined with respect to a given cluster. For finite systems such as molecules it is clear how 

n 

to apply the above formulation [14| . One simply constructs the matrix F from the definition 
in Eq. f|T2|) . where R, R' run over the sites in the molecule. It is not however immediately 
clear how to apply the method to periodic crystals. For this purpose we define a cluster or 
supercell around the unit cell (site) at the origin. A new Wannier orbital centered at the 
origin is constructed as a linear combination of orbitals centered on the sites in the cluster 
as in Eq.Q. For simplicity but without loss of generality let us consider the case of one 
orbital per site (unit cell). First we note that -Fr/ r depends only on the relative distance, 
-^R',R = -^R'-R,o- For R = 0, we obtain the first column F^/ q according to the definition in 
Eq. (|T2l) . We now move to another site R in the cluster and construct a new Wannier orbital 
centered on this site R as a linear combination of orbitals centered on the same cluster sites 
but shifted by R with respect to the cluster centered at the origin. This is to ensure that the 
Wannier orbitals so constructed will be independent of the sites. The cluster centered at R, 
however, has some of its sites outside the original cluster centered at the origin. But there is 
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a correspondence between those sites outside and those sites inside the original cluster: sites 
connected by super lattice vectors are equivalent. This allows us to construct Fr/^r, from 
-Fr',0 in the following way. We search for a superlattice translational vector T such that 



R-o — R.' — R. — T , 



(18) 



is a site in the original cluster centered at the origin. Then the column Fji' n is given by 



A given column of F consists of permuted elements of the other columns. 
B. Results 

Using the procedure described in the previous section we have constructed Wannier or- 
bitals by maximizing U for the 3d transition metals. The cluster consists of the nearest 
and next nearest neighbor sites. As can be seen in Table IIIII the resulting U values are 
remarkably close to the values calculated using the maximally localized Wannier orbitals. 
The results do not change in any significant way when only nearest neighbors are included in 
the cluster. We confirm that the change in U compares favorably with the estimated value 
in Eq. (fT^ . Indeed we have checked that the coefficients of expansion Tn'ii',no is essentially 
unity when H! = 0,n' = n and zero otherwise. This implies that the maximally localized 
Wannier orbitals at the same time to a very good approximation maximize the on-site U 
and form a good basis for the construction of low-energy model Hamiltonians such as the 
Hubbard model. 

To convince ourselves that our procedure is sound, we have performed the following 
calculations. We construct a maximally localized Wannier orbital corresponding to xy orbital 
of the symmetry of SrVOa and calculate U. We also construct Wannier orbitals that are 
deliberately delocalized but span the same Hilbert space as that of the maximally localized 
ones and calculate the corresponding U. A 2x2x2 k mesh is used in the calculations. The 
results are shown in Table HVl As expected, the value of U corresponding to the delocahzed 
orbital is considerably smaller than that corresponding to the maximally localized one. We 
now construct a Wannier orbital by forming linear combinations of both the maximally 
localized orbitals as well as the delocalized ones and maximize U. The orbitals are centered 



-Fr',r — Fro,o • 



(19) 
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TABLE III: The Hubbard U calculated by maximizing the on-site U compared with the values 
obtained using the maximally localized Wannier orbitals. 





uq {<iy ) 


ividA. I e V I 


uu 1 e V ) 




9 AAA^ 1 Q 




o.orj-yjvju 


Ti 








V 


3.216846 


3.216978 


1.3E-004 


Cr 


3.781924 


3.782087 


1.6E-004 


Mn 


3.268736 


3.268788 


5.2E-005 


Fe 


3.619025 


3.619095 


7.0E-005 


Co 


3.097218 


3.097302 


8.3E-005 


Ni 


2.769907 


2.769935 


2.7E-005 



TABLE IV: The Hubbard U calculated by maximizing the on-site U (Max. U) for SrVOa. Starting 
from the maximally localized Wannier orbitals and delocalized orbitals consistently give the same 
maximum value of U . 





Uo (eV) 


Max. U (eV) 


localized 


3.3808733554 


3.3808733554 


delocalized 


3.0292927908 


3.3808733554 



on sites shown in Table IVl The maximum L^'s calculated from the unitary transformation 
of the maximally localized orbitals and the delocalized orbitals consistently agree with each 
other, as they should since the maximally localized orbitals and the delocalized orbitals 
span the same Hilbert space. As in the case of the 3d transition metals, the maximized U is 
essentially identical to the value corresponding to the maximally localized Wannier orbital. 
In Table |V] we show the distribution of weight, |7r^oP, of the orbital that maximizes U 
formed by a linear combination of the delocalized orbitals. The result indicates that the 
original maximally localized Wannier orbital centered at (1 1 1) is the only one that has 
a significant weight at the origin. The corresponding weights for the maximally localized 
orbitals are almost unity when R = and zero otherwise. 

It is remarkable that for the cases considered in the present work the Wannier orbitals 
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TABLE V: The distribution of weight of the Wannier orbitals that maximize U when starting from 
delocaUzed orbitals for SrVOs 



Site 



Weight 



(0 0) 
(0 1) 



4.1E-008 



0.939 



(0 10) 



4.4E-007 



(Oil) 
(10 0) 
(10 1) 
(110) 
(111) 



4.8E-006 



3.1E-007 



6.8E-006 



6.3E-007 



0.061 



constructed by maximizing the on-site U are almost identical with those of maximally 
localized Wannier orbitals. Since the screened interaction Wr is deep around the Wannier 
center, it is reasonable to expect that the maximally localized Wannier orbitals also yield 

a large value of U close to the maximum value. However, the extreme closeness to the 
maximum value is rather unexpected. We have also performed the same calculations by 
maximizing the bare Coulomb interaction and found very similar results. 

V. REAL-SPACE APPROACH TO MAXIMALLY LOCALIZED WANNIER OR- 
BITALS 

Finally we propose a real-space approach of constructing maximally localized Wannier 
orbitals. We construct a unitary transformation on the Wannier orbitals in real space and 
minimize 



We use a combined notation a = (Rn) and the sum is restricted over the sites in a cluster 
or supercell. As in the case of maximizing U consider a small variation 




(20) 



a 



'a ) 



(21) 



18 



S^a = VpTpa, r*p + = , (22) 



a 

= XI [(^"I'^alXa) - {XaKlXaf] 

a 

^^[{Va + S(Pa\rl\iPa + S(Pa) - {(Pa + S(Pa\ra\Va + S(Paf] , (23) 

a 

Ta means that the position is measured with respect to the origin at site a. The change in 
Q to first order in is 

(50 = X [{6ipa\rl\^a) + {^a\rl\SiPa) - 2 {iPa\raM ' {{S^a\raM + {^al^alS^Pa)}] 

a 

= [{^0KMr^a + {'^a\rl\>^f))rf)a-2{ipa\raM ■ {{<^f)\raMTp^+ Mra\^f3)rf3a}] 

a/3 

a/3 

= Yl [{V<^\^1\^P) - i^c^Hl^p)] + 2 [(v?/3|r/3|v:'/3) • (V?a|r;3|(/?^) - MraM ■ {(faKl^Pp)] Tfia ■ 
a/3 

(24) 

As in the case of maximising U we define an anti Hermitian matrix 

Fap = {fa\rl\(Pf3) - (<^a|r||<^/3> + 2 [(^^|r^| ^7/3) • (v?a|r^|v?/3) - MraM ■ {(Pa\ra\^/3)] , 

(25) 

and choose for the steepest descent method 

T = . (26) 

The rest of the procedure is identical to the case of maximizing U. We can write = 
R/3 — Rq + rp implying that fa |¥'/3) — (<^a|i'/3|<^/3)- For a one-band case it follows 
that the third and fourth terms in Fap vanish. The applicability of this scheme depends 
crucially on the feasibility of computing the quantities (falT^alV/s) {fal'^^al^p)- Once 
these quantities are available the minimization process is relatively simple. 
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Following Ref. [5j the spread functional can be split according to 



(27) 







(28) 



The quantity VLi is independent of the unitary transformation. We could equally apply the 
minimization procedure to f2, instead of VL. 



VI. CONCLUDING REMARKS 



We have shown the usefulness of the maximally localized Wannier function as a basis for 
the downfolding procedure. It is found that in transition metals the values of the screened 
Coulomb interaction are in reasonable agreement with those from previous calculations 
based on the LMTO-ASA. The somewhat smaller values of the present calculations may 
be attributed to the more extended nature of the Wannier orbitals compared with the 
more localized truncated partial waves used in taking the matrix elements of the screened 
interaction in the previous calculations. Unexpectedly we have found that for the cases 
we have considered the maximally localized Wannier functions are remarkably close to the 
Wannier function that maximizes the on-site Coulomb interaction. Although we have no 
proof, it is quite likely that this property persists in many other systems. This makes the 
maximally localized Wannier orbitals a very suitable basis for constructing low-energy model 
Hamiltonians. 

We have also proposed a real-space approach for constructing maximally localized Wan- 
nier orbitals, which may be another practical procedure other than the k-space approach. 
The applicability of this scheme, however, remains to be seen. 

There have been many attempts for combining first-principles methods with many-body 
techniques (DMFT and its extensions, path integral renormalization group method [39] 
etc.). The present technique would be useful for the application of these methods to real 
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materials. In particular we have in mind the recently developed GW+DMFT scheme to 
which the present maximally localized Wannier orbitals are now being applied. 
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